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Abstract 

A Gaussian beam method is presented for the analysis of the energy of tlie 
liigh frequency solution to the mixed problem of the scalar wave equation in an 
open and convex subset Q of R", with initial conditions compactly supported 
in il, and Dirichlet or Neumann type boundary condition. The transport of 
the microlocal energy density along the broken bicharacteristic flow at the high 
frequency limit is proved through the use of Wigner measures. Our approach 
consists first in computing explicitly the Wigner measures under an additional 
control of the initial data allowing to approach the solution by a superposition 
of first order Gaussian beams. The results are then generalized to standard 
initial conditions. 
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1 Introduction 

Wc are interested in the high frequency limit of the initial-boundary value problem 
(IBVP) for the wave equation 

n 

Pue = dtUe - 9a., {c^{x)d^^u,) = in [0, T] X n, (1.1a) 

Bu, ^0 in [0,T]xdn, (1.1b) 
Us\t=o = «£, dtUe\t=o = vi in fl, (1-lc) 

where B stands for a Dirichlet or Neumann type boundary operator. 

Above, r > is fixed, £7 is a bounded domain of M" with a C°° boundary and the 
wave propagation velocity c is in C°°{Cl), though this assumption may be relaxed. 
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The initial data depend on a small wavelength parameter e > and we assume that 

ul and arc uniformly bounded w.r.t. e respectively in H^{il) 

and L^{n). (HI) 

We are interested in the description of the behavior of the local energy density 

n 

^\dtUe\^ + c'^\9xjUe\'^ , at the high frequency limit e — > 0, in which case, it is 

well known that this quantity can be computed through the use of Wigner measures. 

The Wigner transform is a phase space distribution introduced by E. Wigner 
[50] in 1932 to study quantum corrections to classical statistical mechanics. In 
the 90's, mathematicians became increasingly interested by these transforms and 
related measures, see for example [29, 33, 34, 35] for the semiclassical limit of 
Schrodinger equations. A general theory for their use in the homogenization of 
energy densities of dispersive equations was laid out by Gerard et al. in [20], see 
also [17, 16]. Wigner measures are also related to the H-measures and microfocal 
defect measures introduced in [49] and [18], see also [6, 1]. Whereas there is no 
notion of scale for the latter measures, Wigner transforms are associated to a small 
parameter tending to zero. In quantum mechanics, this parameter is the rescaled 
Planck constant, while it will be typically the distance between two points of the 
medium's periodic structure for homogenization problems. 

The Wigner transform, at the scale e, is defined for a given sequence (ae,6e) in 
5'(IR")P X 5'(M")P as the tempered distribution 

w,{a„b,){x,0 = (27r)-" / e-'''<a,{x + ^-v)bl{x - ^-v)dv. 

If flg is uniformly bounded w.r.t. e in then ^^[ag] :— Wi;{ai;,a^) converges 

as e goes to in Mp (^S^W^ x W^)j to a positive hermitian matrix measure (modulo 

the extraction of a subsequence), which is called a Wigner measure associated to 
(og) and denoted ?«[&£]. The Wigner measures associated to the solution of the 
wave equation (and hyperbolic problems in general, see e.g. [20, 40]) are related 
to the energy density in the high frequency limit. More precisely, under suitable 
hypotheses (see Proposition 1.7 in [20]), the density of energy associated to the 
solution u'^ of the Cauchy problem for the scalar wave equation converges as £ — >■ 
in the sense of measures to 

where 

1 1 " 

£ («f (i, .)) = -widtu'^it, .)] + w[cdx^u^{t, .)]. 

Moreover, £ (u'^ {t, .)) is the sum of two measures satisfying transport equations of 
Liouville type (see e.g. [20]). 

For the Dirichlet or Neumann initial boundary value problem connected with 
the wave equation, we shall study the same quantity after extending dtUe{t, .) and 
cdxjU^it, .), j = 1, . . . , n, to functions of L'^{K"') by setting dtUe = ladtUe, dj:.Us = 
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^ndxjUir and extending c outside in a smooth way. Hence, we call microlocal 
energy density of Ug the distribution 

1 1 " 

-We \dtUe {t, .)] + We ^cdx/ujt, ■) 

and its high frequency limit the measure 

1 1 " 

£ •)) = 2^ [^tlk(^' •)] + ^ [c9^(i, .) . 

= lj2^e a-nd dxjui = Io^^tjU^ (j = l,...,n) will satisfy the usual assump- 
tions needed in the general context of the study of Wigner measures: their Wigner 
measures are supposed unique and 

vl and dxjul, j = 1, . . . , n, arc e-oscillatory (see (3.23) and (3.24)), (H2) 

the Wigner measures of ^w^^ and (jixjui^ i j = Ij • • • , "-i do not charge the set 
R" X = 0}. (H3) 

Our present study will be restricted to the case where the rays starting from the 
support of the initial data do not face diffraction on the boundary, nor do they glide 
along dft. Therefore, we also assume that 

ul and have supports contained in a fixed compact set of fl independent of e, 

(H4) 

Q is convex with respect to the bicharactcristics of the wave operator, that is 
every ray originating from hits the boundary twice and transversally, 

and 

the boundary has no dead-end trajectories, that is infinite number of successive 
reflections cannot occur in a finite time. 

These geometric hypotheses insure that the only phenomena occurring at the bound- 
ary is the reflection according to the geometrical optics laws. 

Wigner measures for the wave equation in presence of a boundary or an interface 
have been studied by Miller [37] who proved refraction results for sharp interfaces 
and Burq [5] who described their support for a Dirichlet boundary condition. Sim- 
ilar results have been established for other problems [8, 13, 15], in particular the 
eigenfunctions for the Dirichlet problem [52, 19] and for the Neumann and Robin 
problems [7]. All these works are based on pseudo-differential calculus, and in 
particular the use of a tangential pseudo differential calculus. 

In this paper, we present an approach to compute Wigner measures based on 
the Gaussian beam formalism. Therefore, we avoid any use of adapted pseudo- 
differential calculus. Though a Gaussian beam technic requires much more work, 
compared to the above mentioned papers, one advantage is that we are able to give 
asymptotic estimates for remainders terms, which could be useful for numerical 
purposes for instance. 
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Let us recall that Gaussian beams (or the related coherent states) are waves 
with a Gaussian shape at any instant, localized near a single ray [3, 44]. They 
play the role of a basis of fundamental solutions of wave motion and furthermore 
can be used to study general solutions of partial differential equations (PDEs). For 
example, they can help for the understanging of propagation of singularities [44], to 
prove lack of observability [32] and to study scmiclassical measures [41] and trace 
formulas [51, 12]. 

To describe non localized solutions of PDEs, one can use the Gaussian beam 
summation method [24, 10, 25]. The initial field is expanded as a sum of Gaussian 
beams. Each individual beam is computed and the solution is then obtained at an 
observation point by superposing a selection of Gaussian beams. The summation 
strategies are numerous. The sum can be discrete [38, 47, 2] or continuous [30, 31], 
the selection of the beams to be superposed can be done according to several criteria. 
In [4] , a weighted integral of Gaussian beams was designed to build an approximate 
solution of the IB VP (1.1) under an additional assumption (H5) on the initial data 
(see p. 9). See also [27, 28, 43] for recent numerical implementations related to this 
method. 

Gaussian beams seem to be very well suited for the study of Wigner measures. 
Indeed, the Wigner transform of two different beams vanishes when e goes to zero. 
Even better, the Wigner measure of one individual beam is a Dirac mass localized 
on the corresponding bicharacteristic. Thus Gaussian beams act as an orthogo- 
nal family for the Wigner measure. Using these elementary solutions for studying 
Wigner measures is not new, see for example in the whole space domain the work 
by Robinson [45] for the Schrodinger equation, and more recently the paper by 
Castella [9] who used a coherent states approach for the Helmholtz equation. 

As the microlocal energy density of one individual beam is concentrated near 
its associated bicharacteristic, one would expect that the Wigner measure of a 
summation of weighted Gaussian beams will yield easily that the associated weights 
arc transported along the broken bicharacteristic flow (see p. 8 for the construction 
of reflected flows and p. 29 for the definition of the broken flow). Unfortunately this 
result is not immediate as even different beams become infinitely close to each other. 
However, we shall show by elementary computations that this intuition is indeed 
true and that the microlocal energy density of the considered approximate solution 
is transported at the high frequency limit along the broken bicharacterisitic flow. 
Since the asymptotic solution is close to the exact solution u^, we may deduce the 
same consequence for £ {ue{t^ .)). 

The additional hypothesis (H5) consists in assuming that in the frequency space, 
the initial data are supported in a compact that docs not contain (modulo infinitely 
small residues). When studying Wigner measures by the pseudo-differential calculus 
techniques, the frequency behavior of the initial conditions is only controlled by the 
less restrictive hypotheses (112) and (H3). Hence, the assumption (H5) is artificial 
(though not for numerical purposes) and is required only by the Gaussian beam 
summation method we have chosen. However, for e-oscillatory initial data with 
Wigner measures not charging the set R" x = 0}, a truncation of their frequency 
support at infinity and at zero does not affect the energy density of the solution as 
e — >■ 0. By achieving such a truncation, we succeed to derive the transport property 
of the energy density under the traditional hypotheses (H2) and (H3): 

Theorem 1.1 Assume the hypotheses (H1)-(H4) on the initial conditions hold true. 
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Let 8^ = vl ±ic\D\ul. and denote by (pi the broken bicharacteristic flow as- 
sociated to —idt — c\D\ obtained after successive reflections on the boundary 9f2. 



Then 



E {u,{t, .)) = 2 o ^^b'r' + £- o i^lr^) m^x (R"\{0}). 



As mentioned already in our Introduction, this result is well known. But our 
method of proof is able to give more precise estimations than those stated above for 
the Wigner measures. In particular, we have estimations on the Wigner transforms 
of the solutions. 

The rest of the paper is organized as follows. In Section 2, we recall the con- 
struction of first order Gaussian beams and the structure of the asymptotic solutions 
obtained as an infinite sum of such beams. The derivatives of the asymptotic solu- 
tions are then expressed using Gaussian type integrals. We simplify the expression 
of the Wigner transform of such integrals in Section 3, following initial computations 
of [45] in the Sehrodinger ease. We then compute the mierolocal energy density of 
the asymptotic solution by exploiting the expressions of the beams phases and am- 
plitudes and using the dominated convergence theorem. We prove the propagation 
along the broken flow of .) at the high frequency limit, with the help of 

assumptions (112) and (H3) on the initial data. Some complementary results are 
collected in an Appendix, Section 4. 

Let us end this Introduction with a few notations which will be used hereafter. 

A vector x € R"* will be denoted by (.ti, . . . , Xd), the inner product of two vectors 
a, 6 S M'' by a • 6, and the transpose of a matrix A by . If is a subset of M'^, 
we denote E'^ its complementary and \e its characteristic function. For a function 
/ e we let / = Inf. For r > 0, Xr denotes a cut-off function inCg°(R", [0, 1]) 



If no confusion is possible, we shall omit the reference to the lower index x. 
We keep the standard multi- index notations. For a scalar function / € C°°(R!j,C), 
dxf will denote its gradient vector {dxjf)i<j<d and d^f will denote its Hessian 
matrix {dxjdxkf)i<j,k<d- For a vector function g G C°°{R'^, C^), the notation Dg is 
used for its Jacobian matrix {Dg)j^k = dx^gj- If 5 is a function in C°°(Ry x RJJ, C^), 
we denote {Dyg)j^k = dy^gj and {D^g)j^k = dj^kdj- We use the letter C to denote a 
(possible different at each occurence) positive constant. 

For (y^) and (z^) sequences of M+ with e e]0,eo], we use the notation < if 
there exists C > independent of e such that < Cz^ for e small enough. We 
write < e°° or y^ = 0{e°°) if for any s > there exists > s.t. for e small 
enough y^ < CsE'^. 

Finally, if E is in an open subset of M^" and i^^, v'^ are two distributions s.t. 



such that 



Xr(a;) = 1 if \x\ < r/2 and Xri^) = if |x'| > r. 
We use the following definition of the Fourier transform 




lim(i'£ — ly'^) = in i?. 



we shall write 



v'^ in E. 
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2 Tool-box and construction of the asymptotic so- 
lution 



We recall the construction made in [4] of an asymptotic solution as a superposition 
of Gaussian beams and give the expression of its time and spatial derivatives with 
the help of so called Gaussian integrals. 

2.1 First order Gaussian beams 
2.1.1 Beams in the whole space 

Let h+{x,^) = c(a;)|^| and (x*,^*) be a Hamiltonian flow for h+, that is a solution 
of the system 

The curves x^*) of M"+^ are called the rays of P. 

An individual first order (Gaussian) beam for the wave equation associated to 
a ray {t,x*) has the form 

c.,(<,x)=ao(^,x)e^'^(*■-)/^ 

with a complex phase function ^ real- valued on an amplitude function ao 

null outside a neighborhood of (t, x*), and such that 

sup ||Pc.,(t,.)IU^(0) =0(£"), 
te[o,T] 

for some m > 0. 

The construction of such a beam is achieved by making the amplitudes of Pu^ 
vanish on the ray up to fixed suitable orders [44, 23, 32] 

PuJe = (^e-^p{x,dt^,d^ij)ao+e-^t {2dt^dtao - 2c2a,V5xao + ^^^ao) +h.o.t.) e''''/^ 

(2.1) 

where p{x,T,^) — c^(a;)|^p — is the principal symbol of P and h.o.t. denotes 
higher order terms. The first equation is then the cikonal equation 

p{x,dtij{t,x),d,^it,x)) ^0 (2.2) 

on X = x* up to order 2 (see Remark 2.1 in [4] for an explanation of the choice of 
this specific order), which means 

[p (z, dtMt, x), d.^it, x))] U=,* = for \a\ < 2. 
Orders and 1 of the previous equation are fulfilled on the ray by setting 

dti^{t,x') = -h+{x',e) and d^^it,x') = C*. (2.3) 
Choosing ■0(0, x^) as a real quantity, it follows that 

i/'(i,a;*) is real. (2.4) 
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Order 2 of eikonal (2.2) on the ray may be written as a Riccati equation 

+ x')H22{x\e)dli^{t, x') + Hn{x',e) = 0, 

where H ^ ( ) is the Hessian matrix of h^. This nonhnear Riccati 

\^ H21 tl22 J 

equation has a unique global symmetric solution which satisfies the fundamental 
property 

Im cJ^t/i (t, a;*) is positive definite, (2-6) 

given an initial symmetric matrix d'^ij) (O, x*^) with a positive definite imaginary part 
(see the proof of Lemma 2.56 p. 101 in [23]). 

The phase is defined beyond the ray as a polynomial of order 2 w.r.t. (x — x*) 

[48] 

V'(i, x) = ij{t, x*) + ^* . (x - X*) + i(x ™ X*) • dl^{t, x*)(x - X*). (2.7) 

Next, we make the term associated to the power e^^ in the expansion (2.1) vanish 
on (i, X*) 

2dtipdtao - 2c^dxi^dxao + Pipao = on {t, x*), (2.8) 

which leads to a linear ordinary differential equation (ODE) on ao(t, x*). The 
amplitude is then chosen under the form 

ao{t, x) = Xdix - x*)ao(t, x*), 

where d is a positive parameter. The constructed beams are thus defined for all 
{t,x) S ]R"+^ and they satisfy the estimate 

||e-^+iPcj,(t, .)\\L2^n) = 0(Vi) uniformly w.r.t. t e [0,T]. 

Note that Gaussian beams for P associated to the ray (i,x~*) are uje{—t,x). 



2.1.2 Incident and reflected beams in a convex domain 

Assume that c(x) is constant for dist(x, Ct) larger than some constant C > 0. Given 

a point (j/,77) in the phase space T*R", where T*U denotes U x (M"\{0}) if U is 
an open set of R", the Hamiltonian fiow (pQ{y,ri) = {x\^{y,r]),£,l^{y,rj)) satisfying: 

— Xg = c(Xq)|-^, — = — 9a;C(XQ)|^g|, 



is called incident flow. A beam associated to the incident ray (tjXg) is denoted w° 
and called an incident beam. Since we have dependence w.r.t. the initial conditions 
{y,r]), we write the incident beam as 

u;°{t, X, y, 7;) = ao{t, x, y, ,j)e^i^o{t,.,y,n)/e ^ 
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Let TZ be the reflection involution 

7e:T*R"|ao^T*R"|ao 9) 
{X,E)^{X,{Id-2,y{XMXf)E), 

where v denotes the exterior normal field to dfl. We shall only consider initial 

o 

points (y,??) G B = UtgM(^o(^*^) giving rise to rays that enter the domain il at 
some instant. Each associated flow if\^{y,rj) hits the boundary twice. Rcficction of 
^0(^1'?) the exit time t ~ Ti{y,i]) s.t. 



gives birth to the reflected flow (p{{y,ri) = {x\{y , rf) , ^\{y , rj)) defined by the condi- 
tion 

^l^^^^^\y,^l)^no^l^^^'^\y,^). 

Similarly, we also define the refiection time T_i(j/, 77) and the flow Lp^_i{y, rj) by 
reflecting ipQ{y,ri) as follows 



Xq 



\y,v) e an and ±^-'^^'''\y,r,) ■ v (^x];-'^^'''\y,v)) < 0, 

We denote, for fc ±1, the reflected beams by 

iv^it, X, y, r?) = X, ry)e^'/"=(*'--^-'')/^ 

These beams are associated to the refiected bicharacteristics (pj^. Let us introduce, 
for k = 0, ±1, the boundary amplitudes d'L^^^j s.t. 

j=o 

Above, m b denotes the order of B {viib — for Dirichlct and mB = 1 for Neumann) . 
The construction of the reflected phases and amplitudes is achieved by imposing 
that 

1. the time and tangential derivatives of ipk equal at (T^, ccj^'') those of -00 up to 
order 2, 

for fc = ±1. These constraints uniquely determine the reflected phases and ampli- 
tudes, once the incident ones are fixed [44]. If T is sufficiently small, at most one 
reflection occurs in the interval [0, T] and in the interval [— T, 0] for a fixed starting 

o 

position and vector speed {y,rj) £ T*n, and the following boundary estimates are 
satisfied [44] 

||i?(e-t+V(.y,^) + £-^+V(-,y,'7)) IU=([o,T]xaa) = Ole—^-^+i), 
and \\B{e-^+'u;%,y,v)+e-i+'u;-\.,y,Tj))\\H~^^[-T,o]>^on)=0{e-"'--'+i), 

for s > 0. 
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2.2 Gaussian beam summation 

The construction of asymptotic solutions to the IBVP (l.la)-(l.lb) with initial con- 
ditions (1.1c') having a suitable frequency support (see below) is recalled, through 
the Gaussian beam summation introduced in [4]. We focus on a superposition of 
first order beams, for which exact expressions of the phases and amplitudes arc dis- 
played in Subsection 2.2.2. These beams lead to a first order approximate solution, 
close to the exact one up to -^e. Then, the derivatives of the first order solution 
will be approximated by some Gaussian type integrals. 

2.2.1 Construction of the approximate solution 

In [4], we have constructed a family of asymptotic solutions to the IBVP for the 
wave equation for initial data satisfying (HI), (114) and an additional hypothesis 
(115) concerning their FBI transforms. 

Let us recall here that the FBI transform (see [36]) is, for a given scale e, the 
operator : L^{W) L^iM.^'') defined by 

T,{a){y,r]) = c„e-T f a(a;)e"'-(«-^)/^-(«-^)'/(2e)^2;, c„ = a G L^{W'), 

(2.10) 

with adjoint operator given by 

T;(/)(z) = c„e-* / fiy,7^)e"'<--yy'-'---y^"/'-''^dydrj, f e L^R^^'). 

As the Fourier transform, the FBI transform is an isometry, satisfying T*T^ = Id. 
The extra assumption on the initial data needed in [4] is 

l|Tei4llL^(R"xfl^,) - 0{sn and mvi_h^^.^R.) = 0(e°°), (H5) 

where denotes the complementary in W'^ of some ring i?.^ = G vq < \7]\ < 
}, < ro < r 

In general, this assumption may be not satisfied. 

Therefore, we construct a family of initial data ("^ ^.r^^ ' ^e.ro.r^c ) close to 
{ul,v^), satisfying the same assumptions as (HI), (H4) and having FBI transforms 
small in L^(R" x R'^). Letting ro go to and Too go to -l-oo makes these data 
approach [ul^vl) in a sense that will be specified in Section 3.3. In any case, the 
needed convergence is weaker than a convergence since we are interested in the 
study of Wigner measures. 

Let us first truncate T^u^. and T^vl outside i?^ by multiplying them by a cut-off 
lro,r^ G C^(M", [0, 1]) supported in the interior of i?,, 

lro,r^ = Xr„/2(1 " Xin,)- (2.11) 

Lemma 4.5 from the Appendix (Section 4) yields 

rer;7ro,.„c(^)rsi4lU=(K"xfl^) = o{e^), 

and \\T,T:^r,,rAv)Tevi\\L-ij,^^B.^^) = 0(e°°). 
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In order to have data supported in fixed compact sets of fl independent of e, we mul- 
tiply (T*-fro,r^{T^)T,ui,T*-/r„,r^{r,)T,iA by a cut-ofT p g C5"(R", [0, 1]) supported 



in f2, and consider 



It is assumed that p{x) = 1 if dist(2;, suppw^ U suppv^) < C for some C > 0. The 
required estimates 

||T,<,„,,^|U.(K.xfl^) = 0{en and |1T,<,„,,^ = O(e-) (H5') 

are fulfihed since Lemma 4.4 from the Appendix imphes that 

11(1 - p)T:^r„,r^{v)T,ulU. < e-^/- and ||(1 - p)r;7,„,,^ < e'^/'. 

(2.12) 

Using the boundedness of the operator T*^ro,r^T^ from L^(M") to L^(R") and the 
relations 

dyT, = TA^ , d^T; = T*dy^ , (2.13) 

obtained by integrations by parts in the expressions of and T* , one can show 
that the new initial data (u^ ^ ^'^i n, r ) also uniformly bounded w.r.t. e in 

H^n) X 

Let p' be a cut-off of C^(K", [0, 1]) supported in a compact Ky C and satis- 
fying 

p'{y) = 1 if dist(2/, suppp) < C for some C > 0, 

and 7' a cut-off of C^{W, [0,1]) supported in Kr^ C M"\{0} s.t. 7' = 1 on R,,. 
Without loss of generality, we assume that either the incident ray or the reflected 
one propagating in the positive sense is in the interior of the domain at the instant T 
(xj(y, ry) S 51 or x'[{y, rf) £ fi) when y varies in Ky and 77 in R"\{0}. This is always 
possible upon reducing T because the number of reflections for initial position and 
vector speed varying in Ky x (R"\{0}) is uniformly bounded (sec Section 2.3 of 
[4] for similar arguments). And similarly for the instant — T for rays propagating 
in the negative sense. Then, the IBVP (l.la)-(l.lb) with initial conditions (1.1c') 
has a family of approximate solutions u^p^J'^^^ in C°([0, T], H^{n))nC^i[0, T],L^{n)) 
obtained as a summation of first order beams. A general result using a superposition 
of beams of any order was proven in [4], and it reads for first order beams as follows: 

Proposition 1 ([4], Theorem 1.1). Denote for t G [0,T] and x G R" the following 
superposition of Gaussian beams 

= ie-¥+ic„ / p'(y)7'(,y)r,<,^_,^ (y, ij) ( J2 ^e'(.t^ a;, y, rj) 
- X] ^e'i-t^^^y^'nijdydr] 

fe=0,-l 

+ ^£"^+'c„ / p'{yW{rj)e-^T,ul^^^^,^ {y, 77) ( ^ u;^{t, x, y, 77) 

fc=0,l 

+ X! ^ei-t,x,y,ri)jdydr]. 



k=0,-l 
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Above, Wg, are incident Gaussian beams with the same phase ipQ satisfying at 
t = 

^-0(0, X, y,ri)=ri-{x~y)+^-{x- yf (2.14) 

and different amplitudes Og, Og' satisfying 

ao(0, X, y, ?7) = Xd{x - y), [idtijQa^^^ (0, x, y, ?y) Xd{x - y) + 0{\x ~y\). (2.15) 

uif-'^ and Ld^^' denote the associated reflected beams. Then u^^pJ"^^ is asymptotic 
to 

^e.ro.roo cxact solution of the problem (l.la)-(l.lb) with initial conditions 
(1.1c') in the sense that 

sup \\ue.ro,r^ - u^^PP^^^Wm^n) < C {rQ , r^o T) ^/I , 
te[o,T] 

and sup \\dtUe,ro,r^ - dtu'^PPJ^^\\L2(a) < C{ro,r^,n,T) ^/e. 
te[o,T] 

We refer to [4] for further details, and just mention that the proof rehes on the 
use of a family of approximate operators acting from L^(M^") to L^(R"). A simple 
version of the estimate of these operators norms is recalled in Section 4.2 of the 
Appendix. 

2.2.2 Expression of the phases and amplitudes 

In order to compute the first order beams, we begin by analyzing the relationship 
between the incident phase and amplitudes, and the Jacobian matrix of the incident 
flow. A similar relationship involving the reflected phases and amplitudes and the 
reflected flows will be also given. 

The requirement (2.3) for the incident phase implies that 

^ {ipoit, 4)) = dtiJo{t, xl) + dxipo{t, 4) • a;o* = 0. 

Taking into account the initial null value '0g(O,j/) = chosen in (2.14), one gets a 
null phase on the ray 

With the aim of computing ^'^^po{t,XQ), we note that the Jacobian matrix of the 
bicharacteristic Fq ~ D'Pq satisfies the linear ordinary differential system 

iF*, = JH{xieM 
F^ = Id, 

where J = ( 'j-^ ^ j is the standard symplectic matrix. Writing F^ as 



_ I DyXQ -D,,.Tq 

^g — I ^ ^ 



leads to the following ordinary differential system on {Uq, Vq) = {DyX^ + iDj^x^, 

Dyeo + ^Dr,eo) 

jUl = H2iixieoM + H22ixieo)Vl (2.16) 
jV„' = -HnixlQU'o - H,2{xl,eo)Vl (2.17) 
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Note that Fg is a symplectic matrix, i.e. 

{f'.Yjf', = J. 

Using the symmetry of the foUowing matrices 

[DyxlfDyeo,{Dy^)'^D,,eo,Dyxl{D^xl)'^, and i^.^olA^C^)^ 
and the relations 

{DyxlfDr,il - {DyilfD,,xl = Id and Dyxl{D^Q^ - Z?^x* = 
one has 

[Ulf^V^ = (V"o)^f/o, {V^fUl ~ {UlfV^ = 2ild and is invertible. (2.18) 

Putting together (2.16), (2.17) and (2.18) shows that V^o(t^o)"^ is a symmetric 
matrix with a positive definite imaginary part and fulfills the Riccati equation (2.5) 
with initial value ild. Since this is the initial condition for d'^ipo{t, Xq) given in 
(2.14), it follows that 

dlMt,4) = v^iu',)-\ (2.19) 

The incident beams amplitudes are computed as follows. Using (2.3) and the Hamil- 
tonian system satisfied by {xq,£,q), the equation (2.8) at order zero yields the fol- 
lowing transport equation for the value of the amplitude on the ray [23] 

^ (af\t,xl)) + ^Tr (ff2i(4,Co) + ^22(4, ^o)a.Vo(t, 4)) af\t,xl) = 0, 

(2.20) 

which may be written using the matrices Uq and Vq as 
I {af\t,xl)) + ^Tr [{H■2^{xleo)U^ + H^^{xie,)V^) (t/*)-^] at'\t,xl) = 0. 

The time evolution for C/q, see (2.16), combined with the choice of the initial values 
ag(0,y) = 1 and a^^{Q,y) = {-ic{y)\r^\)-'^ from (2.15), yields 

a°(t,4) = (detC/*)-^ and 4 {t,xl) ^ i{c{yM)-^ (dctC/*)^^ . 

Above the square root is defined by continuity in t from 1 at t = 0. 

The expression of the refiected phases i^fc, k = ±1, is similar to the incident 
phase. In fact, since ^ {''Pkit, ^k)) ~ ^ ^nd ipkiTk^x^'') ~ tpQ^Tk, x'^'') because of 
the requirement 1 p. 8, we get 

V'fc(i,4) =0. 

We then apply the general relation between incident and refiected beams phases 
given in Lemma 4.1 from the Appendix, to compute the Hessian matrices of il)±i on 
the rays. As far as we know, the result stated in this Lemma is particularly simple 
enough so that we stated it in the Appendix 4.1. The matrices d'^^±i{t, x^\) can 
also be computed by solving the Riccati equations with the proper values at the 
instants of reflections t = T±i (see eg. [39, 47]). One gets (see Appendix 4.1) 

dlMt,4) = V^iUir^ where C/* = Dyxl + iD^x\ and = DyH + zi?^^!- 
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As ipl is symplectic, {Ul,V^) share the same properties (2.18) as {Uq,Vq) 

{UlfVl = {V^fUl {VlfUl - {UD'^Vl = 2ild and [/* is invertiblc. (2.21) 

The reflected amphtudes evaluated on the rays have an expression similar to the 
incident amplitudes (see Appendix 4.1) 

ag(<,4) = -si (detC/*)"^ and aS'(t,4) = -MyM)-^ (det[/*)"^ for k = ±1, 

where the square root is defined by continuity from i {det Uq''^ ^ at t = Tfe, s = — 1 
for the Dirichlet boundary condition and s = 1 for the Neumann condition. 
We summarize the previous form of the beams in the following result: 

Lemma 2.1 For k = 0, ±1, the incident and reflected beams uj^ have the form 
ujf\t,x)=P,Xd{x-xl)a^^\t)e^^''^', 

with 

/3o - 1, /3i = P-i - -St, 
a,it) = [detUlrK a',{t)^ziciyM)-'[dctUlrK 

= ek ■ - 4) + - 4) • Afc(t)(x - 4), and A,{t) = ~tV^{Ul)-\ 
2.2.3 Gaussian integrals 

It follows that the approximate solution u'^^^^r^ has the form (recall the dependence 
of Gaussian beams w.r.t. variables (y, ry)) 

1 



£ 1 + C 



P'{y)i{v) E Xd{x-xl)hPe,k{t.V,v)e'^'^''^^'''''^''dydri 
ie"T+ic„ / p'ivWiv) Y] Xd{x- x~'^)l3kqe,k{-t,y,rf) 



fc=0,-l 



dydr], 



with 



PeAt^y^il) = ak{t,y,v)£ ^T^ui,,^^,.^{y,r])+a'k{t,y,r])T,v'^,.^^,.^{y,r]), 
and q.^k{t,y,v) = ak{t,y,v)£~^Teui ,,^ ,.^{y,'i]) - a[{t,y,r])Tsvl^^^^^^{y,T]). 

Because of the phases expression given in (2.7), time and spatial derivatives of 
I'iPf^^ may be written as a sum of integrals of the form 

t , ' 1 ' oc 

/ p\yW{v)fe{y,vy - xirr'^^^{t,x,y,T]) 



pii'k{t,x,y,ri)/e 



dydr], j,k = 0,1, |a| < 2, 
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arising from differentiation of .) and U!l{t, .). Other terms of the same form 
originate from derivatives of cj°(— i, .) and i, .). stands for £~^Ti;ul ^.^ 

or T^vl and r^^ arc smooth functions vanishing for jx — a;|,| > d. 
For a function / depending on {t, x, z, 9) e R"+^ x B and fc = 0, ±1, let 

Set K^git) = (filiKy X K^). Let Ilk{t) be a cut-off of Cff (M^", [0, 1]) supported in 
B and satisfying Ilk{t) = 1 on g{t). The volume preserving change of variables 

(.z,0) = ipUy,T]) 
transforms the previous integrals as 

f nfc(Op^''//£-''(a; - zr^fe'^^" ^''^'^''^/'dzde, j, fc = 0, 1, \a\ < 2. 

(2.22) 

We can write the leading terms obtained for j = and a = using Gaussian type 
integrals /^(/i, $) defined as 

for a given phase function $ e C°°(M"^^ x S, C) polynomial of order 2 in x — z and 
satisfying, for t G [0,r] and {z,6) e B 

^{t,z,z,9) is real, dx^{t,z,z,9) = 9, lmdl^{t, z, z,9) is positive definite, (2.23) 

and a given amplitude function h G C°([0, T], L^(R^"g)) supported for every fixed 
t £ [0, T] in a compact of B. By Proposition 3 in the Appendix, one has 

II / h{t, z, 9)x{x - z)e^^'^'-'^-'^-fiy-dzd9Ui < \\h{t, .)||^. 

Noticing that e**/^ is exponentially decreasing for |x — z| > 1, one can use the 
following crude estimate 

II f /i(t,z,6')e^*(*'^'^''')/'^(iz(i6l||i2 < e-'^/'^||/i(i,.)||^2 ^ for a > (2.24) 

J\x — z\>a 

to deduce that Ie{h, $)(t, .) is uniformly bounded w.r.t. e in L^. The same notation 
/e(ft., $) will be also used for a vector valued function h. 

The contribution of the terms (2.22) with j = 1 or |a| > 1 to the derivatives 

of 

'^e'^rj,r^ Order ^/£ as stated in the following Lemma, whose proof is given 
Appendix 4.2 and relies on the approximation operators defined therein. 

Lemma 2.2 (t, .) is uniformly bounded w.r.t. e in i^(K") and satisfies 

dtuZZr^ it^ ^) ^\ Heit, X) - v;-J-t, X)) 

+ 0{y^) m L2(R") uniformly w.r.t. t e [0,T], 



14 



where (v^^) and {v^. ^) are sequences of L?{W^) uniformly bounded w.r.t. e given by 

wt+ = Pkie{-ic{z)\e\iVkff®i'^p;:k^,^^k), 

fc=0,l 
fe=0,-l 

Likewise, c}^;?/"^^''^^ (t, .) is uniformly bounded w.r.t. £ in i^(]R")" and satisfies 

d^uZZr^ it, X) J- X) + V-^'t. ^)) 

+ 0{^/£) m L^(M")" uniformly w.r.t. t e [0,T], 

where (v^ ^) and (v~ ^) are sequences o/L^(M")" uniformly bounded w.r.t. e given 
by 

^Je X! PkleiiO'likP' ® l' P^k'.^k ), 
fe=0,l 

k=0-l 

3 Wigner transforms and measures 

We now compute the scalar measures associated to the sequences (9tu°^^'"^^ (t, .)) 
and (cdxU'^^PJr^it, .))• As \f3k\ ~ 1, the Wigner transform associated to {v^^^{t, .)) 
is a finite sum of terms of the form 

w4Wl,,<fk){t,.), wl^,<i>i)it,.)), 

where k,l = 0, 1, f^^ = cj^ln^p' ® 7' p^'^ and $fc = V'fc ■ As regards the Wigner 
transforms associated to (^cv^ ^{t, .)) , since c is uniformly continuous on R", one has 
by a classical result ([20], p. 8) 

w, {cv+Jt, .),cv+Jt, .)) ^ c^w^v+Jt, .),vXAt, .)) inM2", (3.1) 

and therefore the involved quantities have the form 

cW {Wle. ^k){t, .), WL^ ^l){t, ■)) , 
with fl^ = 9Ukf^'''^k''. 

Similarly, we define for fc = 0, —1 the sequences 5^ ^ = c|0|nfcp' (E) 7' <kJc'^ , which 
are needed when considering the Wigner transform associated with (cw^g(— .)) 
and the cross Wigner transform between (cwj*'g(<, .)) and (ct;^g(— i, .)) , as well as 

, k ^ 

= Ollkp' (81 7' q^k ■ Then, forgetting the powers of s factors, all the previous 
Wigner transforms tested on cut-off functions have the form 

f k ' — ^ ^ 

(3.2) 
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with KgjTg = £ ^ul 
transforms 



and k,l = 0,±1, or after expanding the FBI 




Ks;{w)T^{w')b2''' {z, 6, z', 6' , X, v)e 




■y^(^'^^^^)/^dwdw'dzdedz'd9'dxdv. 



(3.3) 



This type of oscillating integrals is traditionally estimated by the stationary phase 
theorem. For example, this method was suceessfuUy used in [9] for the computation 
of a Wigner measure for smooth data. There the phase was complex and its Hessian 
matrix restricted to the stationary set was assumed to be non-degenerate in the 
normal direction to this set. However, in our case, the amplitude is not smooth as 
no such assumption was made on and v^, and we cannot estimate immediately 
the global integral (3.3) by the same techniques. One possibility of solving this issue 
would be to resort to the stationary phase theorem with a complex phase depending 
on parameters for estimating 



and then study the whole integral involving Ke{w)fe{'w'). 

An alternative method was used in [45], where an integral of the form (3.2) 
associated to the Wigner transform for the Schrodinger equation with a WKB initial 
condition was simplified by elementary computations into an integral over M^". 

Though the method therein faced difficulties in deducing the exact relation be- 
tween the Wigner measure of the solution and of the initial data, we adapt the 
result of [45] to our problem in Section 3.1 and complete the analysis to prove the 
propagation along the flow of the microfocal energy density of u'^^^J'^.^ as e — >■ in 
Section 3.2. The proof is simple and elementary and the computations are made in 
an explicit way. Section 3.3 is devoted to the Wigner measures associated to the 
derivatives of Ue the exact solution of (1.1). 

3.1 Wigner transform for Gaussian integrals 

The sequences (//^g), (/i?.^), {gl.e) ^^i^ (ff^.e) ^"^^ uniformly bounded w.r.t. e in 
L^(M^") and their supports are contained in a fixed compact independent of e. 
Slight modifications of the computations of [45] lead to the following more general 
result: 

Lemma 3.1 Let (/g) and (g^) be sequences uniformly bounded in i^(R^") and 
having their supports contained in a fixed compact independent of e. Let F be an 
open set containing supp/g U supp^g and 5" be phase functions in C°°(R^' x F, C) 
satisfying 



for X e M" and {z,9), {z' ,9') G F, with r<i,,r* G C°°{F,M.) and the matrices 
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H^,Hxji G C°° {F, Mni'C)) having positive definite real parts. Ttien for (f) G C(^{F,R.) 

<We {Ie{fe,^)Je{9e,mA> 

= / (/)(s, cr)/e(s + ^/£r, (T + ^/eS)gl{s - Ver, a - 

A{^, ^)(s, a)e*®-(*^*)("'"''^'*)drdMsda + o(l), 

A($,*)(s,a) ==42*7rt [det (H<i.(s, a) + cr))] 

and 

8e($, *)(s, CT, r, S) —r^{s + y/er, a + y/e5)/e — r,p(s — \/er, a — y/e6)/e 
- 2cr • r/y/e + i{r,6) ■ Q {H^{s,cr),H^{s,cr)) {r,S). 

The matrix Q {Hcp{s, a), H\^{s, a)) and the square root are defined in Lemma 3.2. 

Proof: It consists in two steps. Firstly, the Fourier transform of a Gaussian type 
function is computed explicitly. Then, a Gaussian approximation is used for several 
smooth functions appearing in the Wigner transform integral. 

For simplicity we denote M(a;, z,0) by u and u{x,z',9') by u' when integrating 
w.r.t. z, 6, z', 9' . We also omit the index e in the notation of fe and g^. 
Step 1. Fourier transform. We note that the Wigner transform at point (a;, ^) G 
may be written as 

We{le{f,^),le{g,'f)){x,0 
^^-n^2g-^ f jg*'^ir^/e^ir'^/e+ix-{0-0')/e+i{e'-z'-0-z)/6 



jr ^^-{v+x-z)-His,{v+x-z)/{2e) 

^ ^-{v-x+z'yTu'iv-x+z')/{2e)^ ((2^ _Q_ 

dvdzdz'ddde'. 

The Fourier transform of a Gaussian functions product is given by the following 
Lemma, whose proof is postponed to the end of this Section: 

Lemma 3.2 Let a,b ^ and AL,N G Aidi'C) symmetric matrices with positive 
definite real parts, then 

-p^ (^--(x-a)-M{x-a)l2^-{x-h\N{x-h)l2\^ (-^^ 
= (27r)^ (det(M + 7V))~5 ^-^i-{b^a)l2-{h-aA)-Q{M,N)if,^a,i)IA ^ 

where Q{M,N) is the symmetric symplectic matrix given by 

2M{M + N)-'N ziN - M){M + N)-' 
^{M,i\}~- y i^M + N)-^{N-M) 2(M + 7V)-i 

and the square root is defined as explained in Section 3.4 of [21]. 

Moreover, Q{M,N)A{M,N) = B{M,N) with A{M,N) = .^'^ ) and 

B{M,N) ~ ^ ^ , and Q{M,N) has a positive definite real part 

Re Q(M, N) = 2A{M, N)*-' ^^^^ ^ A{M, NyK 
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Hence 

w, (/e(/,<i>),/e(.9,vI/))(a:,C) 



=42*7r-t£-2" f jg*^ (det(i7$ + H^')) ' ^^^l^-^-^'^l^ 



^i(e+e'-2^)-{z-z')/{2e)+i{e-e')-x/e+i{e'-z'-e-z)/e 

Making the changes of variables 

(z, z') = [s + x/ir, s - V^r), {9, 9') = [g + ^€5,<j - ^e5), 
and writing /+ for /(s + -^/er, cr + ^5) and g_ for g(s — y/er, a — ^5) leads to 

= 42*7r-7e-« /" (dct(i/$+ + iJ* _))"^ ^^r^. + /e-^r^_/e+2^5■(x-s)/^ 

Step 2. Gaussian approximations. Taking the duality product of the Wigner 
transform with a test function (p E C^{F,M.), and after setting {x',^') = {x — s, 
^ — cr)/ y/e, one has 

<w, (4(/,$),/e(g,*))>> 

(3.4) 

Let p'y and be cut-off functions supported in s.t. = 1 on a fixed compact 
containing supp/ and = 1 on a fixed compact containing suppg, and consider 

6e : (x', s, cr, r, (5) t-^ (^(s + ^/ia;', cr + y/e^) 
The r.h.s. of (3.4) may be written as 

<W, (/e(/,<I'),/e(,9,*))>> 

=cl2^n-^ [ (/i(s,a)/+g:l (det(i7<E,+ +i7*_))"^e"-*+/^-*'^*-/^-2«<T.r/y? 

+ c^2^7r-t /" (det(ff<E,+ + i?^_))~V+5le''-*+/"-*'"*-/"-2''^"-/^ 
•^(x' 5')^e(^2(5, 2r, s, cr, r, S)drdSdsda. 

(3.5) 
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Leibnitz formula yields for a multiindex a 

+ P'f+P'9- E C{l3,j)e^d^^,^^,{<j>{s + V~sx',a + V~sa) 



X 



As (s + y/er, a + -v/efJ) varies in supppy and (s — \/er, a — \/e5) varies in supppg, one 
can find by continuity a constant C > s.t. 

Rc i/>i,_) > CId on supp(pj^Pg_). 

Since 

(s,cr) and \/£[r,5) are bounded on supp(p'y^Pg_), (3-6) 
it follows that there exists a constant C" > s.t. 

\d^,^^,h,{x\ e', s, (7, r, (5)1 < Vie-^'("''«')' for all (a;', C', s, a, r, 5), 
which leads to 

\T(^,,^,)be{-25,2r,s,a,r,5)\ < ^/e{l + [r^Sfy-^ for ah {s,a,r,5). 
The second integral in the r.h.s. of (3.5) is then dominated by 

|/+||.g_|(l + {r,5fy~^drd5dsda. 

We deduce by Cauchy-Schwartz inequality w.r.t. s, a that 

<W, (4(/,$),/e(5, *)),</>> 

-42T'7rt j (^(s,ct) (det(i7$+ + i7^_))"V+5le''-*+/'^-*'-*-/"-2ia.r/v/? 



<Ve\\f\\LA\9\\L-, 



where we used AeiQ{Hij,j^, H^, _) = 1 since Q{H^^, _) is symplectic. 

Next, we extend 77$ and i?* outside F as Alf$ + (1 — A)/rf and Ai7* + (1 — X)Id 
by using a cut-off A S C^(M^"', [0, 1]) supported in F s.t. A = 1 on the compact set 
suppp'y U supppjj U supp</), the extended matrices having positive definite real parts. 
The smoothness of these matrices implies by the mean value theorem and (3.6) that 

(det(i7<E,+ + i7^_))"^ - [det {H^{s,a) + H^{s,a))y 
<^/£|(r,(5)| on supp(0/+g:i). 

By symplecticity and symmetry of Q{H^_^_, i?*_), its inverse is —JQ{H,^^, H,^J)J. 
Thus the quantity 

^-(S-r)-Q{H<s. + lU_)-\5-r) _ -{r,S)-Q[H^(s,a);H^{s,a)){r,S) 
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is dominated by 

(r, 5) ■ [Q{H^ + ,H^_) - Q {H^{s, a),H^{s, a)) ] (r, S) 

uG[0,l] 

The positivity of Re and Re Q (iJ$(s, cr), Jy*(s, cr)) and the mean 
value theorem for the matrix function Q (Aiif$ + (1 — X)Id, XH^ + (1 — A)/d) give 
by (3.6) 



-{5,-r)-Q(H^+,H^_)-\5,-r)_ -{r,5)-Q{H^ls,a),H^{s,a})(r,S) 



<^/i|(r,5)|3e-^(-^)^ 



for (s, cr) G suppt/), (s + -v/er, cr + \/eS) S suppp^- and (s — y^r, cr — \/e6) G suppp^. 
It follows that 



< 



V~e\\f\\L49\ 



L2- 



□ 

Proof: [Proof of Lemma 3.2] The matrix M + N has a positive definite real part 
and is thus non-singular. By elementary calculus we have 

(a; - a) • M {x - a) + {x - b) ■ N{x - b) 
= {b - a) ■ M{M + N)-^N{b - a) 

+ (x - {M + N)-\Ma + Nb)) ■ (M + N) {x - {M + N)-\Ma + Nb)) . 

Using the value of the Fourier transform of a Gaussian function (sec Theorem 7.6.1 
of [21]), it follows that 

jr^ ^^-{x-a)-M{x-a)/2^-ix~b)-N{x^b)/2^ 
^-ii-{M+N)-\Ma+Nb)-^-{M+N)-'-i/2 

Writing A/ = 1/2(A/ + TV) + 1/2(M - A^) and TV = 1/2(T\/ + iV) - 1/2(A/ - TV), we 
get the expression with the matrix Q{AI, TV) and the relation 



Q(A/, N)A{M, N) = B{M, N). 



One can easily show that 
B{M,NfjB{M,N) = 



i(Af + TV) 

-i(A/ + TV) 



A{M, NfjA{M, TV), 



from which follows the symplecticity of Q{AI, TV). Then write 

Q{M,N) + Q{M, N) 
=A{M, TV)*"i {A{M, N)*B{M, N) + B{M, TV)* A(A/, TV)) A{M, N)'^ 
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to obtain the value of Re Q{M, N). □ 
From now on, we drop the index e in the notation of v^^, '^^ fte ^tc. for 
simphcity. We fix t G [0,T] and apply Lemma 3.1 with F = B on the sequences 
(ft) Aft) (respectively {fx), {fx)) ^^d the phase functions $fc,<&i for the Wigner 
transforms associated to {yf{t^.)) (respectively (z;+(i, .))). To evaluate the cross 
Wigner transforms between {vf{t,.)) and {v^{~t,.)) (respectively (w^(i, .)) and 
(u~ (— t, .))) , we use this Lemma on the sequences (//""), {g\) (respectively (/^), (ff^))- 

3.2 Wigner measures for superposed Gaussian beams 

We shall prove that the cross Wigner transforms 

We {vt{t,.),Vt{-t,.)) , We (u+(t,.),Wx (-*>•)) 

and 

(4(/L,<I'fc),4(/L,$0) , We (/e(5t'*fc)'^e(gL,*0) 
with k ^ I do not contribute to the microlocal energy density limit £ {u'^ppJ^.^ {t, .)) 

o 

in T*il. We compute 9e($fc, and A($fc, and analyze the transported FBI 
transforms at points {szty/er, azty/eS), which will complete the study of the Wigner 
measures for superposed Gaussian beams. 

Firstly, we note that ||(1 - p' <Si y)UkPe,k\\Ll , = 0{e°°) for fc = 0, ±1. Indeed, 
7' = 1 on Rjj so one gets from (H5') that Tgul , T^v^ have infinitely small 
contributions in L^(R" x supp(l ~ 7')). 

On the other hand, dist (supp(l — p'), suppu^ j.^^ U suppw^ ,,^) > C. Then, 
Lemma 4.3 implies that T^u^ ^^.T^vl^^ have infinitely small contributions in 
L2(supp(l - p') X R"). ThereforcT 

(/,(/^$fe),4(/^$o) 

«A($fc,$0^^ (clalTlk^k") ^ (c|a|n,C/)_ e'^Q^(*-*')drdJ in T*n, 

and a similar relation holds true for [leif^, $fc), Ie{fi, ^i)) ■ 

We start by approaching (c(s)|(t|)_,_ (c(s)|ct|)_ by c(s)^|(Tp in the previous integral 



(3.7) 



and <T+(T*_ by aa* in the integral giving {le{fx, ^k), hifL, 

(4(/^$fc),4(/^,$i)) 

Indeed, these approximations are proved with the help of the following Lemma 

Lemma 3.3 Let (fe), {ge) md $, 4* satisfy the hypotheses of Lemma 3.1. If a and 
/3 are inC^{F,C) then 

W, {Ie{afe, $),/e(/35e,«')) ~ ^fiw, {I,{f,, $),/e(5e, *)) m F. 
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Proof: The proof relies on the use of Taylor's formula on p'y-a and Pg/3, where p'^ 
and Pg are the cut-offs used in the proof of Lemma 3.1 (supported in F and equal 
to 1 on supp/e and supple respectively). □ 
It follows by using (3.7) and (3.8) that 

c^Ttw, (/e(/^-,$,),4(/l,$i)) « (/,(/^ $,.),/,(//, $0) in T*n, 
which leads to 

We {v+{t, .),v+{t, .)) « c^Tr u;, {v^it, .), i;+(t, .)) in T*rt. 

Similarly 

o 

We {vi{~t, .), Vt {~t, .)) w c^Tr We {v~ {~t, .),v^ {~t, .)) in T*f2, 
aT,Awe{vt{t,.),vi{-t,.))^^TYWe{vt{t..).v-{~t,.)) mT*n. (3.9) 

The approximations linking the derivatives of u't'^^^j.^ to v^^ given in Lemma 2.2 
and equation (3.1) lead to 

^We [vt (t, .)] + c^Tr We [v+ {t, .)]+We ^ (-*, •)] + c'Tr We [v- i-t, .)] 
-We {vt{t,.),v^{-t,.)) +C^TrWe {v+{t,.),v^{-t,.)) 
- We {v^i~t, .),v+it, .)) + c'TrWe (v'i-t, .),v+{t, .)) in R^" 

by using the standard estimate (see Proposition 1.1 in [20]) 

I < We{ae,be),(f> > \ < || fle || i2(R„) || 5j || l2(r„) , (3.10) 

for sequences (ae),(6e) in L^(R") and cj) e C^(R^",R). The cross terms between 

o 

v^^ and v^^ cancel in T*Q by using (3.9), leading to 

£ {uZZrJt, ■)) « l^eivtit, •)] + lwe[v-{-t, .)] in T*^. (3.11) 

Thus, we are left with the computation of the Wigner measure associated to {v^), 
computations being similar for (uj ). One has 

We[vt] 

(3.12) 



fc,Z=0,l 



Moreover the inverse of the reflected/incident flow in T*il is a reflected/incident 
flow 

o 

Thus, for {s,(j) G T*il, at most one of the points xZ^.{s,(j) and xZi{s,a) is in il. 
Consequently, the contribution of cross terms between different Gaussian beams in 
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(3.12) vanishes in T*n, and we need to compute only the limits when e goes to zero 
of the following two distributions: 



(3.13) 



Remember that p^^k'' ~ o,k^ £ ^T^u^ 



k --^ k 

+ a'f, T^v^ ,,^ ,.^ , so ^* J. may be written 



as 



--c'{s)\a\^w. 



We 



/e(nfeS^''Te<,„„ ,$fe) 



ic{s)\<j\We { leiUk^k'T.v'e ro , *fc), (HfeS^ Vlr,w| „ ,$ 



(3.14) 



In the remainder of this Section we prove the following Proposition, compute ^* ^ 
and the limit when e of the microlocal energy density of u°-pp^ 



Proposition 2 Let (k^), (t^) be uniformly bounded sequences in L 



2^™"). Then 



[l,iUk^k^T,K,\'i>k),Is{^k^k^T,Te\^k)) ~nlw,{K,,r,)o {^D ' inT*n. 

Above ip\. is extended outside B as the identity. 
Proof: We simplify the integral 



O 



A($fc,$fe) 



Tlk T,T, 



obtained when applying Lemma 3.1 in r*Jl by firstly computing the phase 0^ and 
the amplitude A and then analyzing the transported FBI transforms. Computa- 

tion of 0e($fc,$fc) and A($fc,$fc). We consider (s,cr) G T*i} and start from 



ee($fc, $fc)(s, cr, r, ,5) = -2a ■ r/y/e + i{r, d) ■ Q { Kk{t, s, a),Ak (i, s, a) ] (r, S). 



The particular form of Ak{t) = —iV^{Ul) ^, see Lemma 2.1, induces a similar form 
for the matrix Q Afc (t), A/j (t) 



Q Afc (<),Afe (0 Ffc* 



where Y^* and Z|. are the 2n x 2n matrices 



'^k 



U'k vi 

~k ~h 



and Z{ 



-n n 
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Replacing Ul and by their definitions links and Zl to the Jacobian matrix 



so that 



y; = -iFi J 



-Id Id 

ild ild 



and Zi = J Ft 



k ( -Id Id 

ild ild I ' 



Q{Ak\t),Ak\t) 



As Lp\. o ip_\. ~ Id, one has 



n Fii = Id. 



Combining this relation with the symplecticity of _F^, one gets the following relation 

f~k ~fc ^ 
for the matrix Q I Aj. (t), A^, {t) 



g( Afe-^(t),A, {t)]^{FZlfFZl 



Therefore 



e,($fc, a, r, 5) = -2<j ■ r/^e + i {FZI{s, a)(r, 5)) 



Moving to the amplitude A($fe, <i>fc) — c^2 2 tt2 I det(Afc + A^ ) I , one gets by 



using (2.18) and (2.21) 



Hence 



A($fc,$fc)==c,^^2^"7r2 



detC/* 



Plugging the form of the incident and reflected amplitudes in Lemma 2.1 and using 

'k 



the smoothness of aP on B yields by Lemmas 3.1 and 3.3 



We leiflkak T^Ke ,^k),Ie{^kak T^T^ , $fe) 



.c?,22'Vt 



=:J*_fc(Ke,Te). 



Analysis of the transported FBI transforms. It remains to analyze the most 
difficult terms in the amplitude, which involve transported FBI transforms 



and HfcTeTe ) = (n^ T,Te o ^Zl) [s - V^r, a - ^/^5). 
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Let (j) he a test function in C^(T*fi,M) and -dzl a map of C^(R2"^K2n) ^^la^ 
coincides with (fZk on g{t) U supp^ (see Theorem 1.4.1 of [21]). We use Taylor's 
formula for this map to get for (s ± y^r, a ± y^S) € K^g{t) and (s, a) G supp0 

KDi = ±V^DyXZi r ± V^D,-,xZl S + erf±, 
(Cl) ± = ± ^yCl. r ± V^D.Cl S + erf , 

with 

u)d^dZi{is,a)±uV^{r,S))du, 

u)dyzi ± "V^ '5)) du. 

The change of variables {r',6') = Fzl{s,a){r,S) in i.{k^,t^){s,<7) is thus appro- 
priate. Notice that for (s, cr) e T*n one has the following relations [26] 

DyxZlis^afCti^^'^) - = and -D^x'Ifc (s, cr)^Cfc («, cr) = for w e M. 

In fact, one can show that the derivatives of the previous equations w.r.t. u are 
zero. Besides, the equalities clearly hold true at u = for k = 0, and at m = Tk{s, a) 
for k ~ ±1, as a consequence of (4.9). Hence, it follows that 

<7 - r = dis.a) ■ {DyxZiis,a)r + Dr,xZlis,a)S) = Cl(s, cr) • r', 

o 

which leads in T*n to 

=422"7rt / {nk)+T,K,{xZi + V^r' + er-,+',Cl + V^S' + erl+') 
(n,)_I>7(x:t - Vir' + err',Cl - V^S' + art') 

where 

(rf ,rf)(s,a,r',y) = (r,% r|)(s, a, r, 5). 

In order to use the change of variables (s, cr) = v^K?;, ??) for < J* ^.(Ke, Te), </> >, we 
extend (fil. outside B by the identity and still denote it ipj., making ipl. a one to one 
map from R^" to iflO^'^")- Then UkOipl and (jjoifl belong to C5"(M2",R) and are 



•f (s,a,r,5)= ^ l(r,5)" Al 



l«l=2 
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supported in B. Expanding the FBI transforms gives 

gii7-(2y?r'+£flJ + -ei?;!--2 + 2')/e+f5'-(2l/-2-z'+e-R" + +e-R"")/v^ 

gifl«+-(y+y?r'+e_Rf+-2)-ii?«--(y-v/?r'+ei?,---2') 

g-(!/+v^'-'+eflJ+-2)V(2e)-(a-v/?r'+£i?---2')V(2e) 

where 

(i?^,i?«)(y,r;,r',<5') = (^r,^f)(s,^,r',^'). 
We perform the changes of variables 

(x, u) = , - — -) and y' ^ {y - 

2 e 2 

to obtain 

=c^22"7r^ / + ^u)feix - ^-u)dee''<'-''^ ''dr' dS' dxdudy' dri, 

JR6" 2 2 

where 

=<?!)0V3[,(a; + y/ey\ri){IVkOipl){x + ^/ey' + \/er' + eRl+\r] + ^/eS' + £Rl+') 
(Hfc o(^*fe)(a; + y/ey' - ^/er' + £i?f ?y - ^/e<5' + ei?^"'), 
7e(2;,y',?/,f',<5',w) 

- 7^4-' • (y' -t' + V^Rr' + + ^r'^ + 

+ 7(2;' + r' + - V~£u/2f/2 + - / + x/ii?r' + Viw/2)V2, 

and 

, i?f )(x, y', r?, r', 5') - , + r?, r', 5'). 

Notice that ds{x,y' ,r],r' ,S') converges when £ to 

do{x,T]) = (l)OLpl{x,T])(nkOiplf{x,'q). 

On the other hand, since erj are the remainder terms in the Taylor expansions of 
xZ\{s ± \/er, <7 ± y/eS) at order 2, — is of order ^ and so is — i?^^ , 
leading to 

-(,{x,y',r^y,5',u) ^ 70(2/', r', ,5') = 25' • y' + + 2zr'' + 
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One has 

< J,%(Ks,Te),(/)> 

- ci2^"n^ [ K,{x + ^u)f,ix - ^u)doe'^°e-''T'^dr'd5'dudy'dxdri 

./men 2 2 



< 



(3.15) 

e , . . , e 



\Ke\ix + -u)\t^\{x - -u)dx 
sup \FT^{dee''^' — doe^'"') (x, y', u, r' , 5' , u) \ dr'dS'dudy' . 

X 

Cauchy-Schwartz inequality w.r.t. dx insures that the bracket integral is less than 
II^£||l2||''£||l2- Let us examine the term 

sup I J",, (de 6*'''= — doe'''") (a;, y', u, r\ 5' , u)\dr' d5' dudy' . 

R4„ X 

For fixed y', r', 5' the functions d^ and do are compactly supported w.r.t. {x, J]) so 

sup \T,^[dee'-'^^ - doe'''^°){x,y\u,r',6',u)\ 

X 

< sup I (4e^^= - doe'^") {x, y' , r], /, S' , u) \ . 

Note that |dee*T= - doe'"'°\ is dominated by |4 - do| + Mo| \e''^''-'"'° - l|. The con- 
vergence of c?e when e — > to its limit do is uniform w.r.t. (a;, rj) and so is the 
convergence of 7e to 70 on the support of do. Thus d^e^^' converges to doe'''^" 
uniformly w.r.t. (a;,?/). It follows that 

sup \j^„(d^e''^' — doe^'^°){x,y' ,u,r' ,6' ,u)\ for every y',u,r',6'. 

X e-i-O 

On the other hand, successive integrations by parts give 

'di] = (1 + u^)-" [ L (4e'''=) e-^^-^dT?, 

with L a differential operator w.r.t. 77, of order 2n. Thus, 
sup (dee*'''') (x, y' , u, r' ,S' ,u)\ 



<(l + u2)-»sup max 19" (46*^=) (x, y', r;, r', (5', it) I , 

(a;,r;)l"l<2n 



(3.16) 



for every y',r' ,5',u. The quantities (a; + y/ey',ri) and y/£{r',6') are bounded on 
the support of d^, so i?^*', and their derivatives w.r.t. rj are dominated by 

(r', (5')^. Hence for a given multiindex a, there exists C > s.t. 

1 5^4 1 <C, 

\d^le\ <C\{r',6')\{\{r',S')\ + W + r' + - V~eu/2\ 

+ \y' -r' + ^R^-' - V^u/2\) if |a| > 1, 
for all (a;, y', r;, r', (5') g suppde and u £ M". Thus, there exists C, C" > s.t. 

|9" (4e^'^')| < Ce-C"(y'+r' + y?-R^ + '-\/?u/2)2-C'(y'-r' + v^flr'+V^V2)"-C"r'"-C'5'" 
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for all {x,y' ,r],r' ,6') £ suppd^ and u € R". On the support of d^, ^JeR^^' are 
dominated by |(r',(5')|, which implies for some Co > that 

{2y' + ^eRl+' + Vii?r? > V - Co|(r', 
Hence, if \y'\ > Co|(r', e-c-'^^y'+v^^^+'+v^^r')' < e^'^"?'". Otherwise, 

g_CV'^-C'5'= < g-C"s/' = -C"r'^-C"5'^^ jj^ g^jj ^^^^g^g^ ^-^^^^ ^^^ists C", C" > S.t. 

for every x,y' ,r],r' ,S' ,u and e e]0,eo] with some Eq > 0. Using this in (3.16) leads 
to 

sup \j',j{d,e'-'-){x, y', u, r', S' , u)\ < {1 + u2)-"e-^^"-^''"-^''", 

X 

and repeating the same arguments for sup ((ioe''''°)| gives 

X 

sup \ j-rj{di:e''^' — doe^'^") (x, y', u, r , 6' , u) \ < {I + u^Y "e 

X 

for every y',u,r',S' and e €]0,eo]- By the dominated convergence theorem, one 
obtains 

sup\J^^{d,e^^' - doe'''"){x,y',u,r',S',u)\dy'dudr'dS' 0. 

X £— ^0 

From the inequality (3.15) concerning the distribution f^{Ki,,Ti;), one finally has 
by plugging the expressions of do and 70 

n f £ S 

< Jl,k{^e, Te), >= c^22"7r 2 / K^{x + -u)f^{x - -u) 

e^'^'-y'-y"-^''"-^"e-''^-"dr'd6'dxdudy'dr] + o(l). 
Integration w.r.t. r' , S' , y', ry yields 

< >/lfc('«£,T-e),0 > 

= (271")"" / J^n{nlo(fl(j)Oipl){x,u)Keix + ^u)fs{x - ^u)dxdu + o{l). 

The integral in the r.h.s. is exactly the Wigner transform of (keiT^) tested on 

nloipl(j}Oipl. □ 

We are now able to compute the measure /i* f. given in (3.14) by using the 
previous Proposition and the Lemma 4.7 

/i*,,. « nl {w, - tc\D\ui^^^^^J) o {^i)-' in T*^l. 

Recalling the relation between the Wigner measure and the FBI transform (see 
Proposition 1.4 of [19]) 

/ \TM^\^edydi^ ^ <w[a:\,e (3.17) 
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for G C5"(R^",R) and (oe) a uniformly bounded sequence in L^(R"), it folfows 
that We [vi^ro,r^ ~ ^'^\-^\'^i,ro,r^] ~ ^ i-^y ^ '^j?)'^ o^' equivalently 

w;e [<,„,,^ - tc\D\ui,^^^,J o i^iy' « in (At,(t))^ 

Since 11^ = 1 on g{t), one deduces 

By summing over fc = 0, 1 and letting e — >■ 0, we get 

fc=0,l 

For u e [— T, T] and (y, 77) G A'j, x (R"\{0}), the incident and reflected flows are 
related to the broken bicharacteristic flow associated to —idt — c\D\ as follows: 

{'P-i{y,v) if w < r_i(?/,77), 
Voiy.v) if r_i(y,r/) < w < ri(y,ry), 
(Pi{y,v) iiu>Ti{y,7]). 

We extend at times of reflections arbitrary. We define (p^ in {fl\Ky) x (M"\{0}) 
by successively reflecting the rays at the boundary. As only one incident/reflected 
ray can be in the interior of the domain at a fixed time t G [^T, T] 

</)o^*= ^ 00^*, ini^, xR"\{0}. 

A:=0,1 

It follows that 

— 1 ° 

•)] = W K,ro,r^ " ^^iDWi^ro^r^] » (fl) m T*n. 

The computations for are similar. One has just to replace the index k = 1 
by fc = — 1 and by in (3.13) and to repeat the same techniques. If we 

denote '^f^ro,r^ = "^Lro.r^ ^ , then one gets 

— 1 ° 

w{vt{-t, •)] = w o (</J^*) in r*r2. 

Using these results in (3.11) as e — s- leads to 

(3.18) 

3.3 Proof of the main Theorem 

A consequence of the estimate (3.10) is 

I <w{ae,he),9 > I < limsup ||ae||L2(o)limsup ||fee||i2(Q), (3.19) 
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for (oe), (be) uniformly bounded sequences in L'^{W) and 9 e C^{T*n,R). Apply- 
ing this estimate to the difference between the derivatives of the exact and approxi- 
mate solutions of the IB VP (l.la)-(l.lb) with initial conditions (1.1c'), one deduces 

the measures associated to (^dtUe^ro,rSj ^^^d ^9xMe,r-o,r^^ and gets by (3.18) 



o {ipl) ^ in 



T*n. 



Remark 1 Gaussian beam summation of first order beams allows to compute the 
microlocal energy density of the solution of the IBVP (1.1) as e — > 0, under the 
hypotheses (H1),(H4) and (H5) on initial conditions. Summation of higher order 
beams may imply asymptotic formulas for the Wigner transforms and thus for the 
energy density. Higher order terms in the expansion of the Wigner transform were 
studied for instance in [14] and [42] for WKB initial data. 



Let us now study the microlocal energy density for the problem (1.1) when e — > 0, 



by making the data {uI j.^ j.^,vIj.^ j.^) approach {ul,v^). The contribution of the 
sets {t] G K", \r]\ > roo/4} and {77 e M", \r]\ < 4r-J^ where jro,r^ ^ 1 (remember 
the definition of 'yro,r^ in (2.11)) to T^ul,T^vl is controlled asymptotically by the 
assumptions (H2) and (H3). 

Set Tf =vl± ic\D\ul and denote 0* = (/)o</j*. Then (/)* e C^(M2",M) and one has 



<£(u,(i,.)),0>-- <«;[T+] 



> 



<- |< w \T+ 



-I- |< It; [dfujyt, .)] - w 

n 

cdxjU^it, .) 



>l + 2 l<^fco-oo] -^[T^r] >l 



> 



< w 



> 



> 



< w 



[T- 



(3.20) 



Wc use (3.10) to get 



<limsup||T+,,^, 
<limsup \\vl ~ I 



1l2(r,.) limsup (||T+ 



limsup llwg 

6-i-O 



,rJ\\Hi{n)- 



Similarly, by (3.19) 



< W \dtuAt, .)] - W dtUe^ro,r^ jt, .) 



> 



<limsup \\dtUs{t, .) - dtUe^ro,r^{t, OlU^Cn) 



limsup ||afUe(t, .)|1l2(J2) + \imSVLp\\dtUs^ro.r^{t,-)\\L^n) , 
e-)-0 e~>0 
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and for j = 1 , . . . , n 



< w 



> 



<linisup Wdx^u^it, .) - dxjUs.ro,r^ (t, ■)\\mn} 

e-i-O 

( limsup \\dxjU^(t, + limsup WdxjU^.ro.r^ [t, ■)\\L^{n) 

The solution of the IBVP for the wave equation is given by a continuous unitary 
evolution group on the space H^{fl,dx) x L^{^l,dx). Hence 

\\dtU,{t, .) - dtU,^ro.r^{t, .)\\L2^n) < -«e,ro,r^lU2(0) + II " l^^ro 1 1 ^ ^ (O) , 

\\dx^U,{t,.) ~ dx^U,^ro,r^{t,-)\\mn) < hi - vlro.rJlL^n) + 1 1 ^ W^,^^, || i (n) , 

for j ~ 1, . . . ,n. Finally, by using (3.18), the estimate (3.20) is simplified into 

I < £(u,(t,.)) ,0 > ~l<w [T+] ,0-* > _i < u; [TJ] I 

2 2 (-3 21) 

<limsup||i;^ - vI^^ ,,J\l2(^^^ +limsup||u^ - ul^^^^Jlmin}- 

e-i-O ' ' e-i-0 

We therefore need to estimate the difference between initial data (1.1c) and (1.1c'). 
We start by the initial speed. By the exponential decrease of T*^ro,r^T^vl on the 
support of 1 — p (see (2.12)), one has 

Because T* is bounded on L2(M2") ^ L'^(W) and T^T, = Id 

\\v^-T*'^Tevl\\L2{Rr.) < 11(1 - Xroo/2)rewi_||L2(H2,.) + \\Xr^ /2X4roTevl\\L^R^r.) . 

® ® 

Firstly, Lemma 4.2 yields 

{®f = ||c„(2^)-t£-t (1-x.^/2(^)) I -Fz4(Oe'«-^-^''-^«^'/''^'rfeili2(K-). 
It follows by Parseval equality that 

{®f = cle-'-i I (1 - Xr^/2{ri)f\:Fvi{0\\-^^--i)'/^-didri 

+cle-^ [ {l-Xr^/2m'\^vlm'e-(^-^i^'/^d^dr^. 
The first integral in the r.h.s. is exponentially decreasing, which leads to 



limsup ® < limsup ( / \J'vl{£_)\^d^ 



\"'kCI>' 

Secondly, as dist (suppw^ , supp(l — p)) > 0, one gets ||(1 — p)Tevl\\L2^jg.2n-^ < e~'^l^ 
by Lemma 4.3 and thus 

limSUp( ®f limSUp||/?(j/)x,.^/2(?/)X4ro(?7)rewf ||i2(R2„). 

E-J-0 E-S-0 
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It results from the relation (3.17) applied with = vi that 



Because w is a regular measure, assumption (H3) yields 

Va > 0,3/o(a) > s.t. w Iv^] ({|CI < ^o(a)}) < a. 
One deduces, for 4ro < lo{a), that 



lim sup ® < \/a, 

e-i-O 



which leads to 



limsup||u^ - vl^^^ ,.^\\L2f^n-) < limsup / + Va- 

e^O £^0 \J\e^\>r^/8 — J 

For the analysis of — it^ H^{fl), we begin by estimating the spatial 

derivatives of the difference. It follows by using the relation (2.13) when differenti- 
ating the inverse FBI transform that 

The term involving the derivative of p is exponentially decreasing by Lemma 4.4. 
Since the FBI transform of a derivative is the derivative of the FBI transform by 
(2.13), one has to estimate — pT*^ro,r^Ti,dxjUl\\L2{^[-iy Employing the same 

previous techniques yields for j = 1, . . . , n 

limsup||9:,^.u^-9:r,Ue,ro,r^||L2(f2) < limsup / \^ {dx,ul{m\'^dn +v^, 

if 4ro < and w dx ul ({|^| < lj{a)}) < a. Set vq — \ min lj{a), then the 



Poincarc inequality yields the same bound for limsup||u^ — , 

e-yO 

Coming back to (3.21) we deduce that 

<£iu,it,.)),^>-^<w [T+] ,</>-*> -1 < z„ [T7] 



> 



lim sup 



l?l>'-o=/8 



(3.22) 



The assumption (H2) of e— oscillation means by definition that 
limsup / iJ'U^m'de^^ 0, 



lim sup 

e-i-O 



£|«I>-R 



for j = 1, . . . , n. 



(3.23) 
(3.24) 
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Since the l.h.s. of the estimate (3.22) does not depend on a nor Too, one deduces 
by taking the hmits a — >■ and r^a — > oo that 

£ {u,it, .)) = Iw [T+] o i^-y' + [TJ] o ((^*)-' in T*n. 



4 Appendix 

4.1 Reflected first order and higher order beams 
4.1.1 Higher order beams 

Higher order beams, possibly with more than one amphtudc, can be constructed to 
satisfy better interior and boundary estimates. In this case, the eikonal equation 
(2.2) must be satisfied up to order i? > 2 on the rays. If r > 3, the equations 

{p{x, dti^, a,^)) (i, X*) = 0, |a| = r, (4.1) 

give systems of linear ODEs of order 1 on (5"'0(t, x*))|^|^^ with second members 
involving lower order spatial derivatives of the phase. In fact, the key observation 
is the equality 

drp{v')dtd2il^{t,x')+di_p{^') ■ dxd^^{t,x') 
^2c{x')\e\dtd:^{t,x') + 2c2(a;*)e* • d^d^^it, x*) 

=2c(x*)|C*||(a,"^(t,x*)), 

used for \a\ = r to eliminate the r + 1-th order derivatives of ■0 in equation (4.1). 
To summarize, the requirements 

dMt,x') = ~c{x')\e\, d,ij{t,x') = e, 

p {x, dt4'{t, x), dxipit, x)) — on X = X* up to order R, 

uniquely determine the spatial derivatives of '0 on the ray up to the order R under 
the knowledge of their initial values on (0,a:°). We refer to [44] for further details. 



4.1.2 A general relation between incident and reflected beams phases 

By (2.19), the Hessian matrix of the incident beam's phase is related to the Jacobian 
matrix of the incident flow. One can prove that its higher order derivatives are also 
related to the higher order derivatives of the incident flow. Computations exhibiting 
such relations can be found for instance in the Appendix of [39]. We shall give a 
nice relation between an incident phase -01110 and the associated reflected phase -^ref 
for beams of any order. This relation is intuitive true on geometrical grounds and 
it provides with the derivatives of the reflected phase up to order R, which might 
be useful in applications of Gaussian beams. 

Consider the following auxiliary function linking ip\ to ip^ for any flxcd time t 

{x,0^ vf'^''''^ ono^l^^^-'\x,0- 
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For a given point (x,^) £ B, si{x,^) is its "image by the mirror" dH.. For in- 
stance, Chazarain used this type of auxiliary functions in [11] to show propagation 
of regularity for wave type equations in a convex domain. 

By the Implicit functions theorem, Ti is C°° on the open set B and so is Si. 
Since o si satisfies the same Hamiltonian equations as (p\ and ip^^^'^'^\x,£,) = 
ifQ^^'^'^^ o Si{x,^) for {x,^) e S, one has 

ip\ =flo Si. 

Besides, noticing that Ti{(p^) = Ti — t, one has also 

= si o (4.2) 

(fQ and (p\ are symplectic C°° diffcomorphisms from B to B [22], and so is si. One 
can define a similar auxiliary function s_i : B ^ B s.t. Lp*_i = tp^ o s_i and 
= s_i o (/Sg for t e R. 
Let us introduce the components of si as 

si = (r. A). 

For m e N, /, g functions in C°° (w^ x (M^\{0}), , uo e K" a fixed point and 
V e C°°{R'^,C^) a phase function s.t. V{uo) e M^\{0}, we introduce the notation 

f{u,V{u)) " g{u,V{u)), 

U — Uq 

to denote that the formal partial derivatives of / (u, V{u)) and g (u, V{u)) up to the 
order m coincide on uq. The differentiation here is viewed formally, since V may 
be complex valued out of uq, which makes f{u,V{u)) and g{u,V{u)) not defined 
for u Uq. However, on the exact point uq, one can always use the formula of 
composite functions derivatives to get a formal expression of the derivatives. We 
will use the same notation 

f{t,x,V{t,x)) ^ git,x,V{t,x)), 

X — X^ 

for functions /, g e C°° x x (]R^\{0}), C^) and phase function V € C°°(Rt x 
R^,C^) s.t. for t € R, V{t,x*) £ R^MO} to denote that the formal partial deriva- 
tives of / (t, X, V{t, x)) and g it, x, V{t, x)) w.r.t. x up to order m coincide on (t, x*) 
for all t G R. We will be sloppy with respect to the notation of the dependence of 
the phase V on its variables. 

Consider an integer R > 2 and an incident phase V'inc satisfying 

R 

dt^inc{t,xl) = -C{xl)\^l\, dj:^inc{t,xl) = ^-^d p{x, dt^inc, d^i^inc) ^ 0. 

As a particular case, the phase tpo is obtained by setting R ~ 2 and choosing its 
initial value on the ray as zero and its initial Hessian matrix on the ray as ild. 

Let ipj-ci G C°°{M.t X R", C) be the reflected phase associated to V'inc, that is the 
phase satisfying 

dtiprci{t,x\) = -c(x*)|^*|, dxil}^af{t,x\) = £,{ and pix^dtAcfyd^ipref) ^ ^ 0, 

x—x\ 
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and having the same time and tangential derivatives as i/'inc at the instant and the 
point of reflection (Ti, Xq^) up to the order R. 

Since (p^ and the reflection TZ conserve c(x)|f| (see (2.9)), one has for every 
{x, C) e S and T e R* 

P {r{x, r, A(x, ^)) = p{x, r, 

Thus 

which implies, by construction of ipinc 

R 

p{r{x, do,1pinc),dt'^inc,Mx,da:i}inc)) t 0. (4.3) 

Compare this with the equation 

R 

p {r{x, 92;-0inc), 9t-0rof (t, r{x, 9:r^inc)) , d^lprcf it, r{x, d^lpinc))) ^ t ^ 

resulting from the construction of iprei and (4.2). This suggests the following Lemma 
Lemma 4.1 

dtll^ref{t,r{x,dxli)^nc)) StV'mc and da:^ref (t, r{x , doo-ip^nc)) ^^\{x , dxll}i„c) ■ 

A similar result linking the reflected phase associated to the ray (i, .tI* ) to -^inc can 
be established. 

Proof: The strategy of the proof is the following: we consider a phase function 
9 satisfying the relations announced in Lemma 4.1 and we prove that 9 fulfills the 
eikonal equation on the reflected ray up to order R and has the correct derivatives 
at the instant and point of reflection. This proves that 9 coincides with the reflected 
phase on the reflected ray up to the order R. 

Denote r (x, dxipinc{t, x)) by Q{t, x) or simply by g if no confusion arises and let 
us first verify that for a fixed k>l there exists a phase function 9 e C°°(Mt x R", C) 
s.t. 

8.^,0) ^^^^\{x,d,t.nc)- (4.4) 

Let A{t,x,0 = D^r{x,£) + D^r{x,£)dli:;,,^{t,x) and B{t,x,0 = D^X{x,^) + 
D^X{x,^)dlil;inc{t,x). Then DxQ{t,x) = A{t,x,da;ipinc), D^[X{x,dxiJinc{t,x))] = 
Bit, X, a^Vinc) and for v € C°°(Rt x R^, C^) one has 

Dx {v{t, g)) X ^ D^v{t, g)A{t, x, S^^V'inc)- 
Hence, 9 exists if A(t,Xo,^o) ^-On singular and 

S(t,a;,5a;V'inc)^(i,a;, ^j^Vinc)"^ 'x^^ (A(t, S^j^'inc)^) ^ B(t,X,dxTpir,cf ■ (4.5) 

From (4.2) one gets 

A{t, Xq , ^q) {DyXQ + iD^jXo) = DyX^ + iDj^x^. 
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Since Lp\ is symplectic, the matrix ^ ^^^t ^ symplectic. This imphes in 

particular the relation 

D.,ei{Dyx\f ~ Dy£X{D,,x\f ^ Id 

and the symmetry of Dyx\{Dnx\)'^ . Thus, ker(£',,a:;^)^ n ker(Z?j,a:^)^ = {0} and at 
the same time, 

{Dyx\ + iD^x\){Dyx\ + iD^x\)* = Dyx\{Dyx{f + D^x\{Dr,x\f. 

This proves that Dyx\ + iDnx\ is invertible and so is ^(t,Xo,^o)- On the other 
hand, 

.4 \ / D^r D^r \( Id 
B )-[ D^X )[ dlt^n, 

Since JM = DsJjDsi, the symplecticity of si leads to 

M^JM = J. 

Hence 

and the requirement (4.5) is fulfilled. 

The relation (4.4) fixes the derivatives of dtdxO on {t,x\) up to order k — 1. 
Indeed, using the compatibility condition 



f {t,X,dxtpinc{t,x)) I 



d 
dt 

=dt [f {t, X, dxipincit, X))] l^^^t + dx, [f {t, X, dxll^incit, x))] \ 



on the maps (t,x,^) i-> dxd {t,r{x,(,)), {x,£_) i-> X{x,£_) and their derivatives yields 
recursively by (4.4) 

dt [dx9{t,g)] 'x^^ D^X{x,dxipinc)dtdxipinc- 

Thus 

dtdxO{t,g) + dl6{t, g)D^r{x,dxilJinc)dtdx^inc ""^^ D^X{x, dx^inc)dtdxipinc- 



k— 1 

Using the relations d'^0{t, g) x ^ (iJA^^) [tjXjdxtpinc) and (4.5) in the previous 
equation yields 

dtdxO{t, g) 



fc-i 



D^Xix^dxTpinc) - ((^^) ^-B^) {t,X,dx'lpinc)D(^r{x,dxll^inc) 



dtdx'ipu 
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Since 

A-D,X -B-D,r=^ ^^^^^ ) " M^JM ( ) = 

it follows that A{t,x,dx4'inc)'^ dtdxQ{i, q) dtdx4>inc- Note that 

dx {u{t, g)) ^ Ait, X, d^^incfdMt, e) for u € C°°(]Rt x M^, C), (4.6) 
so one gets 

dx {dtO{t, g)) 'x^^ dtdxipinc- 
Setting dt9{t,x\) = dtipinc{t, Xq) implies then that 

dt9{t, g) ^ dt^inc (4.7) 
Putting together (4.3), (4.4) and (4.7) shows that the phase 6 satisfies 

p{g,dte{t,g),dx9{t,g)) 

under the further assumption k > R. 

Let 7r(i, x) = p {x, dtO{t, x), dxd{t, x)). Since dx (7r(t, g)) {t, Xq) ~ and 
A{t, Xq,^q) is non singular, it follows by (4.6) that dxTT{t, x{) is zero. More generally, 
for TO > 1, the formula of composite functions' high derivatives yields 

n n 

^^.1 [7r(i,e(t,a;))](t,x^) ^ ^^.1 ■■■dx,^^it,x{)Y[Aj^,^{t,xl,Q 

fc=i 

where depends on derivatives of tt on (i,x*) of order lower than m. For 

m < R, the l.h.s. is zero so one can show recursively on < R that 9f 7r(t, x* ) = 0. 
One thus has the following eikonal equation on 9 

p{x,dt9,dx9) X 0. 

x—x^ 

To compare the time and tangential derivatives of 9 and V'inc at (TijxJ^^), let us 
introduce a C°° parametrization of a neighborhood U of .t|^^ in 9r2 

where A/" is an open subset of M"^^, (T{J\f) = and tr is a diffcomorphism from J\f 
to U. For .T G M" close to a;^\ we may write a; = cr{v) + VnV {<7{v)) , with v ^ Af 
and Vn G K. Denote <y{vi) = x^^ and set 9h{t,v) — 9{t,a{v)) and (V'inc)f, (i> = 
V'inc (i, ""(i")) the phases at the boundary near x'^^ . Since r{X, "E) = X for (X, S) G 

T*M"|an, it follows that 

(t,«) = (Ti,«i) 
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which implies by (4.7) that 

k 

dtOb ^ dt (^inc)t • 

(t,«) = (Ti,«i) 

Similarly X{X, S) = S - 2 (S • i^{X)) i^{X) for (X, 5) e r*R"|ao, leading to 

Da(wf A(a(z;),5,Vmc(i,a({)))) x Da(i))^a,Vinc (t, . (4.8) 

(t,t)) = (Ti,{ii) 

Since di,9b{t, v) = Da{v)'^dx6 [t, (t{v)) and a similar relation holds true for di, (^inc)^, 

k 

one gets from (4.4) and (4.8) that dijOb ^ dy (V'inc)h • Hence 6b and {'il'inc)b 

(t,fi) = (Ti,f)i) 

have the same time and tangential derivatives at (Ti , vi ) from the order 1 to the 
order k + I. 

If we assmne that 9{Ti, Xq'^) = ■(/'incCTi, Xq^), then 

0b ''^l . (^inc)fc, 

and 6 satisfies all the requirements that determine the reflected phase associated to 
tpinc and concentrated on (t, x\). The phases 6 and V'rcf ai'c thus equal on {t, x\) up 
to the order R. □ 

4.1.3 First order reflected beams' phases and amplitudes 

Lemma 4.1 gives at order one 

dlMt. xl) {DMxieo) + D^r{xlil)dlMt, 4)) 
^D^X^xlil) + D^X{xl iDdlMt. A)- 

One obtains by plugging the expression (2.19) of d^tpo{t,XQ) 

dlMt,x\) {D.,r{xl,eo)U'o + D^r{xl,eo)V,') = i?,A(x* , e*)f/* + i?eA(.T* , e^)Fo*. 

From (4.2), it follows that 

dlMt,4) = V^iUl)-' where C7* = Dyxl + zD.xl and = Dy^ + zD.CL 

and a similar relation holds true for d'^ip-i(t,x'Li). 

The reflected amplitudes evaluated on the associated rays satisfy transport equa- 
tions which are similar to (2.20) and may be written as 

^^(^af\t,xi))+lTi-[{HM,a)Ul + H,,{xlQV^) {Ul)-'] af\t, xl) ^ 0. 

One can obtain a similar equation to (2.16) on Uf. involving H2i{x\t^\) and 
-^22(2:^, Cfc) I by using the relation ip\ = ip^ o Sk- On the whole 



k 
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In order to find the relationship between U'^'' and Uq^ for fc = ±1, we differen- 
tiate the equality x^'' — x^'' 



where the square root is obtained by continuity from 1 at t = T^;. 
On the other hand, for k = ±1 

+d'l„,^ ^b{x,d^^o)a'^ + b{x,d^i)k)ao, 

where b denotes the principal symbol of B. Thus, the condition 2 p. 8 required for the 

construction of the reflected amplitudes implies that Cq (T^ , a;Q ) = sajj (Tfc , ccq ) , 
with s = — 1 for Dirichlet condition and s = 1 for Neumann condition. 

k ^Tfc 

and compute the derivatives of Tk from the condition Xq^ G 90 

dy.,Tk = ^ {Dy^.x^^fi^ix^") 

i^xo^" -lyix^")) 

to get after elementary computations 

U^" - [id - 2K4'= H^" f) U^" ■ (4.9) 

Hence 

4{t,xi) = -si {dctUiy^ and 4'{t,xl) = sic{yMr' (det[/*)"^ for k = ±1, 



w 



here the square root is defined by continuity from i[dct Uq''] 2 at t = T; 



4.2 Approximation operators 

We briefly recall a simple version of the integral operators with complex phases 
used in [4] and the estimates established therein. We then use these results to prove 
Lemma 2.2. 

For t € [0,r], let Kz,e{t) be a compact of M^" and consider the set 

El = {{t, X, z, 0) G [0, T] X R3", (z, 0) e K,^g{t), \x-z\< 1}, 

which we assume compact. Let $ be a phase function smooth on an open set 
containing Ei and satisfying (2.23) for t € [0,r] and {z,9) G Kz,e{t). Then there 
exists r[$] e]0, 1] s.t. 

Im$(i,x,z,6') > C{x ~ zf for t e [0,T], {z,e) e Kz^g{t) and |a; - z| < r[<I>]. 

Let G C°°([0,r] X R3«^c) satisfying 

for t e [0,r],L(t,x,z,6') = if (2,6*) ^ K,^g{t) or |x - z| > r[$], 

e^S^Je is uniformly bounded in L°°{[0,T]x R^") for every 1 < j < n and /s e N. 

(4.10) 

If O" {le{t, .), $(t, .)/e) denotes, for a given multiindex a and t G [0,T], the operator 
[O" mt,.),^t,.)/e) h] {x) 

h{z, e)L{t, X, z, d){x ~ z)"e**(*'^'^'^)/"dzd6', h e L2(k2„)^ 
then, under the previous hypotheses on $ and Z^, we have the following estimate: 
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Proposition 3 ([4], Lemma 3.3) 

||0"(;e(t,.),*(i,-)/e)IU^(K-)-.L^(K") < uniformly w.r.t. i e [0,T]. 

This estimate allows to prove Lemma 2.2. 

Proof: [Proof of Lemma 2.2] Consider the integrals (2.22) giving the derivatives 

~k 

of u^PPJ^^ and fix j, k and a. The transported phase ipk is smooth and satisfies 
by (2.3), (2.4) and (2.6) the properties (2.23) for t € [0,T] and {z,9) e K^,g{t). 

We fix some €]0,1] so that Im^k'' {t, x, z,e) > C{x - zf for t G [0,T], 

{z,B)&K\,(t) and \x^z\<r^k\- 

For t G [0, T], Efcp' (g) 7' {t,z,d){r^- ^ {t,x,z,9) depends smoothly on its variables 

~k 

and vanishes for ja; — z| > d or {z,9) ^ K^g{t). Hence, upon ehoosing d < r[-0fc ], 
the amplitude H^p' ® 7' \r^_a) satisfies the properties formulated in (4.10). Let us 



- k 



k k 



check if Ig/e ~ IjsTevl ,.^ ,.^ ,1b£~^T^uI j.^^ ,.^ is uniformly bounded in L^(M^"). 
Clearly T^v^ is, and the property holds true for e~^T^ul by Lemma 4.6. 
One can then use the approximation operators O" to write the integral (2.22) as 

£-T+J / Hfep^'*' f.^ri^Jix - z)°'e''^'''''^dzde 

=£-T+JO" (^n,p^''(^ (t,.),V^'(t,.)/e^ 1b/,'. 

The estimate established in Proposition 3 yields 

||£-T+J [ Hfep^^'" /e'(^\x - zre'^^'l'dzdeWL^ < e^+^. 

Hence, only (rpg) contributes to dt^xUe^ror^-. the residue being of order y/e. One 
has 

ro,o{'t,x,y,r]) = ^Cnl3kdt,a:ipk{t,xl)xdix - xl)a[\t,y,r]), 

and by (2.3) 

9tVfc(t,4) = -c{xiMi dxMt,xi)=ek- 

It follows that 

=0^""^^" / J2 i-i)Pkciz)\0\xdix - z)Ukit,z,e)p' ®y {t,z,e) 



^''(t,z,6')e*'^'= ^''^'''^^/"dzde 

£ 4 r„ / 

2 

fc=0,-l 

qTk^i-t, z, e)e''^'''~^-*''='''''')/')dzd0 



1 3„ /• --^ fe 

-s ^ Cn > i/3fcc(z)|6l|xd(a;-z)Hfe(-t,z,6l)p'(8)7' (-t,z,( 

•^K'" fe=0,-l 
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in L2(k«)^ uniformly for t e [0,T] and 

1 3 f ^ 

+ -e-Tc„ / V ip,kOoXd{x- z))lu{'~t,z,B)p' ®i {-t,z,e) 



fc=0,-l 



0(x/?) 



in uniformly w.r.t. i e [0,T]. 

One can get rid of the cut-off Xdix — z) appearing in dtu'^J'^^^ {t, x) 
and dc,^u'^PPJ,.^{t,x) by using the estimate (2.24). ' ' °° □ 

4.3 Results related to the FBI and the Wigner transforms 
Lemma 4.2 For u in L^(M") 



Proof: The equality is proven by Parseval formula. □ 

Lemma 4.3 ([4], Lemma 2.4) Let a he a positive real, E a measurable subset ofR" 
and K C M" a compact set s.t. dist(i4r, £') > a. If u £ L^(R2) is supported in K 
then 



\\TM\lhe.^^j = c„e-T|ll^(y)^,(x)e-(--^)'/(2^)|U. ^ < e-'^'/^^-) ||«|| 



Proof: The proof consists of writing the FBI transform as the Fourier transform 
w.r.t. X of some auxiliary function and using Parseval equality. □ 

Lemma 4.4 Let 9 he a cut-off o/C^(K^,M), E a measurable subset of and 
K C M" a compact set s.t. dist(i4r, £') > 0. If u E L'^(]R") is supported in K then 

\\T:e{ij)TM\LHE)<e-'"'\\u\\L^^„y 

Proof: The kernel of lET*e{ij)T^lji L^{R^) ^ L^{Rl) is 

ke{w,x) = e^^cllEWKix) I 0(,y)e'''-(---)/^-(y--)V(2e)-C^-a)V(2e)^y^^ 



^ lB(w)lif(x)e^"(27r)-"J-0(^^)e-("^"')'/(4"). 
For w e M", one has by Cauchy- Schwartz inequality 

\ke{w,x)\dx < l|^(?|U2(R„)(27r)-"e-t (^j^ lE{w)lK{x)e-^-^^" /^^^'Ux^ 
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Similarly, \k^{'w,x)\dw is dominated by e so one gets by Scliur's Lemma 

\\T:9iv)T,u\\mE^,) < e-c/^||M|U.(Kj). 

□ 

Lemma 4.5 Let E be a measurable subset o/ M" and K C M" a compact set s.t. 
A\si{K,E) > 0. If e IS a cut-ojf of C(^{K;^,R) supported m K then 

Proof: Consider the operator : ^^(IR^;^ ^^(M^"^) defined by 

It is easy to compute its kernel 

e, y, 7?) = e-"l£;(C)l^f (,7)e^(«+''H--y)/(2^)-(-y)V(4e)-(?-^)V(4e)_ 

Hence, \he{x,^,y,r])\dxd^ < e"*^/^ and J^^^ \he{x,(,,y,r])\dydri < e^'^/". For 
u € L^(M") , it follows by Schur's Lemma that 

ll^erew||i2(R2„ ) = \\T,T;9{rj)T,u\\L2(^s,^y,E^) < e-^/^||reii||L2(R2„j 

<e-^/^h|U2(R„). 

□ 

Lemma 4.6 ([4], Lemma 3.4) \\e-^T,ul ,.^^^J\L2(^^2n-^ < 1. 
Proof: Differentiating (2.10) w.r.t. yj, < j < n, yields 

£^dy^{Teui,ro.rJ = T-Vj^^^T.ul ,.^ ,.^ - c„£"T / (x)e" ^ (y^ - Xj) 

The l.h.s. is bounded in because dy-{T^u[_^^^^^) = Tg{dxju[_^^^^^). The second 
term of the r.h.s. is the Fourier transform of a bounded function in i^, thus it can 
be estimated using Parseval equality. One gets 



Thus lle^ST^jTjW^^^ 11^2 < 1 and consequently by (H5') 

Hence roolU^ ^ V^- Reproducing the same arguments on the equality 

dy,{Teui^ro.rJ T-V]£^^T^ui.ra,r^ " CnS'"^ / (e^^ W^^^,,^^^) {x)e~^ {y^ - Xj) 

^iv{y-x)/e-iy-xf/{2e)^^ 

leads to Ijuf |li2(R„) <e. □ 
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Lemma 4.7 Let (og) and (be) be two sequences uniformly bounded in i^(R") and 
H^{W') respectively. If £~^bi: is uniformly bounded in L^(R"), then 

Weia,, \D\be) « \^\we{ae,e"\) on R" x (M"\{0}). 

Proof: Let (/> be a test function in Cg°(R" x (M"\{0}),R) and denote = \D\be. 
We use another expression of < Weiae, Ce),(j) > exhibiting the Fourier transform of 

Ce- 



< We (fle, Ce ),(/)>= (27r) " / T^(l){x — -V , v)ae{x)ce{x — ev)dvdx. 
Since J-^cj) is rapidly decreasing 



sup 



By Cauchy-Schwartz inequahty w.r.t. dx 

I (T^(f){x - - u)) ae{x)ce{x - ev)\dvdx < £\\ae\\L^2\\ce\\L^. 

It follows that 

< We{ae,Ce),4> >^ {27r)~" / (j){x, S^)e~™'^ ae{x)ce{x - ev)dvdxd^ + o{l) . 



Integrating w.r.t. v leads to 

<Weiae,Ce),^>=i2Tr)-''e-'' 



{x, Oe-"-«/"ae(x) J"ce(^/£)dxdf + o(l), 



and replacing Fce{£,/e) by e ^|^|^6£(C/£) ends the proof. 



□ 
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